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Abstract 

We present a measurement of the rate of gluon splitting into bottom quarks, g —>■ bb, in 
hadronic decays collected by SLD from 1996 to 1998. The analysis was performed by 
looking for secondary bottom production in 4-jet events of any primary flavor. A topological 
vertex mass technique was used to tag the two jets with the smallest angle between them as 
b/b. We obtained a rate oi g bb per hadronic event to be (3.07±0.71(stat.) ±0.66(syst.)) x 
10~^ (preliminary). 
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1 Introduction 



The process of the sphtting of a gluon into a heavy-quark pair is one of the elementary 
processes in QCD but is poorly known, both theoretically and experimentally. 

The rate is defined as the fraction of hadronic events in which a gluon splits into a 
bb pair, e'^e' qqg —>■ qqbb. The value of g^i is an infrared finite quantity, because the 
6-quark mass provides a natural cutoff, hence it can be safely computed in the framework 
of perturbative QCD However the rate is sensitive to the as parameter and to the 
6-quark mass, which results in a substantial theoretical uncertainty in the calculation of gi,i. 
The limited accuracy of the (7^5 prediction is one of the main sources of uncertainty in the 
measurement of the partial decay width = r(Z° 66)/r(Z° — > qq) [0, In addition, 
about 50% of the B hadrons produced at the Tevatron are due to the gluon splitting process, 
and a larger fraction is expected to contribute at the LHC A better knowledge of this process 
can improve theoretical predictions of heavy-flavor production at such colliders. 

This measurement is difficult experimentally. The cross section of g ^ bb is very small 
even at Z° energies, since the gluon must have sufficient mass to produce the bottom-quark 
pair. There are huge backgrounds from Z° — > bb whose magnitude is about a hundred 
times larger than the qqg qqbb process. Moreover the B hadrons from g —>■ bb 

have relatively low energy and short flight distance and are more difficult to distinguish 
using standard vertexing. So far, the only three measurements of (7^5 have been reported, by 
DELPHI and ALEPH [|. 

Here we present a new measurement of gi,i based on a 400k Z^-decay data sample taken in 
1996-98 at the Stanford Linear Collider (SLC), with the SLC Large Detector (SLD). In this 
period, decays were collected with an upgraded vertex detector, wider acceptance and 
better impact parameter resolution, thus improving considerably the 6-tagging performance. 

2 The SLD Detector 

A full description of the SLD and its performance have been described in detail elsewhere 
Only the details most relevant to this analysis are mentioned here. 
SLD is well-suited for the measurement of g ^ bb due to two unique features. The first 
is that the SLC, the only linear collider in the world, provides a very small and stable beam 
spot. The SLC interaction point was reconstructed from tracks in sets of approximately 
thirty sequential hadronic Z^ decays with an uncertainty of only 5//m transverse to the 
beam axis and 32/im (for bb events) along the beam axis. Second is the upgraded vertex 
detector (VXDS) [Q, a pixel-based CCD vertex detector. VXDS consists of S layers with 
SOOM pixels and each layer is only O.S6% of a radiation length thick. The measured r0 {rz) 
track impact-parameter resolution approaches 11/im (2S/im) for high momentum tracks, 
while multiple scattering contributions are 40/im/(px sin^^^ ^) in both projections {z is the 
coordinate parallel to the beam axis and p± is the momentum in GeV/c perpendicular to 
the beam line). With these features, topological vertex finding gives excellent 6-tagging 
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efficiency and purity. In particular, the efficiency is good even at low B-meson energies, 
which is especially important for detecting g bb. 

3 Flavor Tagging 

Topologically reconstructed secondary vertices [0 are used by many analyses at the SLD for 
heavy-quark tagging. To reconstruct the secondary vertices, the space points where track 
density functions overlap are found in 3-dimensions. Only the vertices that are significantly 
displaced from the primary vertex (PV) are considered to be possible B- or D-hadron decay 
vertices. The mass of the secondary vertex is calculated using the tracks that are associated 
with the vertex. Since the heavy-hadron decays are frequently accompanied by neutral 
particles, the reconstructed mass is corrected to account for this fact. By using kinematic 
information from the vertex flight path and the momentum sum of the tracks associated 
with the secondary vertex, we calculate the Pr-corrected mass Mp^, by adding a minimum 
amount of missing momentum to the invariant mass, as follows: 

VTX 

Here Myrx is the invariant mass of the tracks associated with the reconstructed secondary 
vertex and Pt is the transverse momentum of the charged tracks with respect to the B-flight 
direction. In this correction, vertexing resolution as well as the PV resolution are crucial. Due 
to the small and stable interaction point at the SLC and the excellent vertexing resolution 
from the SLD CCD Vertex detector, this technique has so far only been successfully applied 
at the SLD. 



4 Monte Carlo and data Samples 

The measurement uses 400k events collected from 1996 to 1998 with the requirement that 
the VXDS was fully operational. 

For the purpose of estimating the efficiency and purity of the g —>■ bb selection proce- 
dure, we made use of a detailed Monte-Carlo simulation of the detector. The JETSET 7.4 
Q event generator was used, with parameter values tuned to hadronic e~^e~ annihilation 
data [^, combined with a simulation of B hadron decays tuned to T(45') data [|l^ and a 
simulation of the SLD based on GEANT 3.21 |]ll|. Inclusive distributions of single-particle 
and event-topology observables in hadronic events were found to be well described by the 
simulations ||T2[. Uncertainties in the simulation were taken into account in the systematic 
errors (Section 0). 

Monte-Carlo events are reweighted to take into account current estimates for gluon 



splitting into heavy-quark pairs 0, 0. The JETSET at SLD predicts gi,i = 0.14% and 
gcc = 1.36%, and we reweighted them so that g^^i = 0.273% and gcc = 2.58%. A Monte-Carlo 
production of about 1200k Z ^ qq events, 1000k Z ^ bb events and 480k Z ^ cc events 
are used in order to better evaluate the efficiencies. 
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Besides the signal events, hereafter called B, two categories of background events exist: 

• Events which do not contain any gluon splitting into heavy flavor at all, hereafter called 
Q events; and 

• Events in which a gluon splits to a charm quark pair, named C events. 

5 Event Selection 

The two B hadrons coming from the gluon tend to be produced in a particular topological 
configuration, which allows one to discriminate the signal from background. We select g ^ bb 
events as follows: 

• Require 4 jets in the events; 

• Require b tags in two jets selected in a particular configuration; and 

• Apply additional topological selections to improve the signal/background ratio. 



Jets are formed with energy- flow particles, using the Durham jet-finding algorithm [14 
with ycut = 0.008, chosen to minimize the statistical error. The 4-jet fractions for the B, 
C and Q events predicted by the simulation are about 32%, 18% and 5.3%, respectively. 
The overall 4-jet rate in the data is (5.976 ± 0.044)%, where the error is statistical only. 
In the Monte-Carlo simulation the rate is (5.678 ± 0.002 ± 0.068)% where the first error is 
statistical and the second is due to the uncertainty in the simulation of heavy-quark-hadron 
physics. The two jets forming the smallest angle in the event are considered as candidates 
for originating from the gluon splitting process g bb. The selected jets are labeled as jet 
1 and jet 2, where jet 1 is more energetic than jet 2. The other two jets in the event are 
labeled as jets 3 and 4, where jet 3 is more energetic than jet 4. 

Jets containing B-hadron decay products are then searched for by making use of the 
information coming from the vertex detector, using the topological vertex method. We 
require jet 1 and jet 2 to each have a secondary vertex. No tag is applied to jet 3 and jet 
4. After topological vertexing, about 300 events are selected. The selection efficiency for 
g ^ bb is expected from Monte Carlo to be 6.6% while the signal/background ratio is 1/5. 
67% of the background comes from Z ^ bb events, 21% from g ^ cc events and remaining 
12% from Z qq {q ^ b) events. 

In order to improve the signal/background ratio, we use topological information. Firstly, 
many bb background have one 6-jet splitting into 2 jets so that the two found vertices are 
from different decay products from the same B decay. The two vertex axes tend to be 
collinear. Figure |I| shows the angular distribution between vertex axes in jet 1 and jet 2. 
Half of the bb background peaks at cos ^12 ~ 1. In order to remove bb events, we require 
-0.2 < COS012 < 0.96. 

Secondly, the variable | cosq;i234|, where 01234 is the angle between the plane II12 formed 
by jets 1 and 2 and the plane 1134 by jets 3 and 4, is used to suppress the bb background. 
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Figure 1: Angular distribution between vertex axes in jet 1 and jet 2 (0.9 < 0056*12). Points 
indicate data, open box signal, hatched boxes are backgrounds. 



Figure ^ shows the distribution of | cos 0:1234 1. This variable is similar to the Bengtsson- 
Zerwas angle |T^, and is useful to separate g ^ bb events because the radiated virtual gluon 
in the process qqg is polarized in the plane of the three-parton event, and this is 

reflected in its subsequent splitting, by strongly favoring g qq emission out of this plane. 
Events with | cos 01234! > 0.8 are rejected. 




COSCX,234 

Figure 2: The distribution of the cosine of angle between the plane ni2 formed by jets 1 and 
2 and the plane 1134 formed by jets 3 and 4, for data (points) and Monte Carlo (histogram). We 
reject | cos 01234] > 0.8 



Thirdly, the b jets coming from a gluon tend to have lower energy than the other two jets 
in the event. We require the jet-energy sum of jet 1 and jet 2 to be smaller than 36 GeV. 

Finally, c jets have lower P^^-corrected mass than b jets. Figure § shows the greater of 
the Pr-corrected mass determined for jet 1 and jet 2 after the above cuts. Many g —>■ cc 
events are in below 2.0 GeV. Hence we require maximum P^-corrected mass to be greater 
than 2.0 GeV to remove g ^ cc events. 



6 Result 

After requiring all the above mentioned cuts, 62 events are selected in the data. Background 
events are estimated to be 27.6 using Monte Carlo, where 63% of the background comes 
from Z bb events, 27% from g ^ cc events and the remaining 10% from Z ^ qq {q ^ b) 
events. Table |l| shows the tagging efficiencies for the three categories of events, where the 
errors are statistical only. From these efficiencies and the fraction of events selected in the 
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Figure 3: Maximum Pj^-corrected mass distribution between jet 1 and jet 2 after jet-energy-sum 
cut. Points indicate data, open box signal, hatched boxes are backgrounds. 



Source Efficiency (%) 



B 
C 

Q 



3.86 ±0.52 
0.10 ±0.02 
0.73 ±0.05 



Table 1: Efficiencies after all cuts for the three categories. Errors are statistical only. 

data fd = (2.14 ± 0.27) x 10~^, the value of (7^5 can be extracted as: 

_ /d - (1 - gcc)eQ - gccec 

9bb — • 

- eg 

The measured value of the gluon splitting rate into bb pairs is 

= (3.07 ± 0.71) X 10-3, 

where the error is statistical only. 



(1) 



(2) 



7 Systematic Error 

The efficiencies for the three event categories are evaluated by Monte-Carlo simulation. The 
limitations of the simulation in estimating these efficiencies lead to an uncertainty on the 
result. The error due to the limited Monte-Carlo statistics in the efficiency evaluation is 
Aqi^I = ±0.44 X 10"^. This uncertainty comes mainly from the efficiency to tag Q events. 

A large fraction of events remaining after the selection cuts contain b and c hadrons. 
The uncertainty in the knowledge of the physical processes in the simulation of heavy-flavor 
production and decays constitutes a source of systematic error. All the physical simulation 



6 



parameters are varied within their allowed experimental ranges. In particular, the b and c 
hadron lifetimes as well as production rates are varied, following the latest recommendations 
of the LEP Heavy Flavour Working Group |16|. The uncertainties are summarized in Table 

i 

The simulation of the signal events is based on the JETSET parton shower Monte Carlo, 
which is in good agreement with the theoretical predictions In order to estimate the 
uncertainty on this assumption, we have produced 10,000 g ^ bb events using GRC4F [|r^ 
at the generator level. The signal tagging efficiency mainly depends on the energy of the 
gluon splitting into bb. This efficiency function, computed with JETSET, is reweighted by 
the ratio of GRC4F to JETSET initial distributions to obtain the average efficiency. A 
systematic error of ±0.09 x 10"'^ is estimated from the difference in ee from the two Monte- 
Carlo models. 

The dependence of the B efficiency on the 6-quark mass has also been investigated at 
the generator level. Events are generated using the GRC4F Monte Carlo, which is based on 
a matrix element calculation including 6-quark masses. The variation of the B efficiency is 
computed as done for JETSET, using the GRC4F spectrum for 6-quark masses from 4.7 and 
5.3 GeV/c^. The uncertainty is estimated to be 0.06 x 10~^. 

The uncertainty in the ratio of the g —>■ cc background events, Ag^c = ±0.40%, gives the 
error Ag^i = ±0.09 x lO'^. 

There is about 5% discrepancy of 4-jet rate between data and Monte Carlo in our i/cut- 
The uncertainty due to the discrepancy is estimated by increasing background events in the 
Monte Carlo to be Ag^i = ±0.14 x 10"^. 

Charged Monte-Carlo tracks used by the topological vertex tag are smeared and tossed to 
better reproduce distribution of data. Uncertainties in the efficiencies due to this smearing 
and tossing are assessed by evaluating the Monte-Carlo efficiencies without the smearing and 
tossing algorithm. The difference in the g^^i result is taken as systematic error. The errors on 
gf,l due to the tracking resolution and efficiency are then estimated as Agi,i = ±0.26 x 10~^ 
and = ±0.29 x 10^^, respectively. 

Table |^ summarizes the different sources of systematic error on gi,i, and the total sys- 
tematic error is estimated to be 0.66 x 10~^. 



8 Summary 

A measurement of the gluon splitting rate to a bb pair in hadronic Z° decays collected by 
SLD has been presented. Excellent SLC and VXDS performance provides advantages not 
only for 6-tag efficiency but also for topological selections. The result is 

g^b = (3.07±0.71(stat.) ±0.66(syst.)) x 10"^(preliminary). 

where the first error is statistical and the second includes all systematic effects. 
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Source 


^9bb ( 


10-^) 


Monte Carlo statistics 


±0, 


.44 


b hadron lifetimes 


±0, 


.01 


b hadron production 


±0, 


.07 


b hadron fragmentation 


±0, 


.12 


b hadron charged multiplicities 


±0, 


.11 


c hadron lifetimes 


±0, 


.01 


c hadron production 


±0, 


.03 


c hadron charged multiplicities 


±0, 


.08 


Energy distribution oi g ^ bb 


±0, 


.08 


b quark mass 


±0, 


.06 


Qcc 


±0, 


.09 


4-jet rate discrepancy 


±0, 


.14 


IP resolution 


±0, 


.09 


Track resolution 


±0, 


.26 


Tracking efficiency 


±0, 


.29 


Total (Preliminary) 


±0, 


.66 



Table 2: Systematic uncertainties on g^^g. 
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